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ABSTRACT 

The magnetosphere of a rotating pulsar naturally develops a current sheet 
beyond the light cylinder (LC). Magnetic reconnection in this current sheet in- 
evitably dissipates a nontrivial fraction of the pulsar spin-down power within a 
few LC radii. We develop a basic physical picture of reconnection in this environ- 
ment and discuss its implications for the observed pulsed gamma-ray emission. 
We argue that reconnection proceeds in the plasmoid-dominated regime, via an 
hierarchical chain of multiple secondary islands/flux ropes. The inter-plasmoid 
reconnection layers are subject to strong synchrotron cooling, leading to signifi- 
cant plasma compression. Using the conditions of pressure balance across these 
current layers, the balance between the heating by magnetic energy dissipation 
and synchrotron cooling, and Ampere's law, we obtain simple estimates for key 
parameters of the layers — temperature, density, and layer thickness. In the co- 
moving frame of the relativistic pulsar wind just outside of the equatorial current 
sheet, these basic parameters are uniquely determined by the strength of the re- 
connecting upstream magnetic field. For the case of the Crab pulsar, we find them 
to be of order 10 GeV, 10 13 cm -3 , and 10 cm, respectively. After accounting for 
the bulk Doppler boosting due to the pulsar wind, the synchrotron and inverse- 
Compton emission from the reconnecting current sheet can explain the observed 



1 Center for Integrated Plasma Studies, Physics Department, University of Colorado, UCB 390, Boulder, 
CO 80309-0390 

2 Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 



-2 - 



pulsed high-energy (GeV) and VHE (~ 100 GeV) radiation, respectively. Also, 
we suggest that the rapid relative motions of the secondary plasmoids in the 
hierarchical chain may contribute to the production of the pulsar radio emission. 



Subject headings: Gamma rays: stars — Magnetic reconnection — Radiation 
mechanisms: non-thermal — Relativistic processes — pulsars: general — pulsars: 
individual (Crab) 



1. Introduction 



Simple theoretical reasoning shows that a plasma-filled magnetosphere of a rotating magnetic 
dipole, e.g., an axisymmetric pulsar, cannot remain closed beyond the ro tator's Light Cylin- 
der ( LC) radius, -Rlc = cf2 -1 , where Q is the pulsar spin frequency (IGoldreich fc Julian 
19691 ). Instead, the magnetic field lines extending beyond the LC have to open up, in- 



evitably forming an equatorial current sheet (CS), see Fig. [TJ The current sheet splits into 
two separatrix surfaces at the so-called Y-point, which is probably located at the LC. This 
magnetospheric configuration has been studied extensively by a variety of numerical ap- 



proaches, including relativ i stic fo rce-free Grad-Shafranov equation (jContopoulos et al.lll999 



Gruzinov 



20051: 



Timokhin 



tions ( Komissarovl 20061: McKinnev 



20061). time-dependent force-free electrodynamic (FFE) simula 



2006 



Spitkoyskv 



2006 



Kalapotharakos fc Contopoulos 



b: 



20091 : Kalapotharakos et al 



2012 



Li et al. 



2012 



Petri 



tivistic magnetohyrodynamic (MHD) simulations (IKomissarov 



2012 



Parfrev et al 



2006 



2012). and rela- 



Bucciantini et al.ll2006l ) 



and is now well understood. Moreover, its non-axisymmetric generalization, co rresponding 



to an oblique pulsar, has also been obtained by time-dependent FFE simulations ( Spitkovskv 



2006 



2012a). 



Kalapotharakos &: Contopoulosll2009l : Kalapotharakos et al.ll2012l ; iLi et al. 



2012: 



Petri 



However, despite this remarkable progress in recent years, we cannot regard the above steady 
state solutions as completely realistic. In particular, all of these numerical models suffer from 
one common problem — they cannot properly treat the singular equatorial CS and have to 
resort to various artificial procedures to keep the codes stable. While the exact choices for 
these assumptions and procedures do not seem to have drastic effects on the overall gross 
structure of the magnetosphere, including such basic global parameters as the open poloidal 
magnetic flux fraction and the total spin-down power, they do affect strongly the behavior of 
the equatorial CS, especially near the Y-point. We would like to stress that what happens in 
this region is important because it likely involves dissipation of a fraction of the outpouring 
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electromagnetic energy via magnetic reconnection (see below). This is important for two 
reasons. First, reconnection alters magnetic field topology and t hus affects the structure of 
themagnetosphere and of the pulsar wind (IContopoulosI 120051 ; IContopoulos fc Spitkovsky 



20061 ); in particular, the bursty, non-stationary character of the reconnection process marked 



by p lasmoid eject i on is likely to lead to quasi-periodic motions in the magnetosphere (see, 



e.g 



Contopoulosl 120051 ; ISpitkovskyl l2006t iBucciantini et al.l 120061 ). Secondly, reconnection 
in the equatorial CS (or "ballerina's skirt" CS in the case of an oblique pulsar) means 
that a sizable fraction of the pulsar spin-down power is dissipated locally near the LC, 
instead of contributing to the Poynting flux powering the pulsar wind at larger distances. 
A significant fraction of this dissipated energy is very likely to be radiate d away promptly 



see § [2D and thu s powers the observed high-energy gamma-ray emission (|Lvubarskii 



20001; 



2012a 



Aronsll201ll.l2012l; iBai k, Spitkovskyll2010l ; IContopoulos k, Kalapotharakod 



2010: 



1996 



Petri 



Arka fc Dubus! 20121 ). Figuring our the character (e.g., the spectrum) of this emission 



cannot be done completely within the framework of relativistic MHD or force-free models 
but requires a closer look at the basic plasma physics inside the reconnecting and radiating 
equatorial CS. 



Determining these physical conditions is the main goal of the present paper. This includes 
both the description of the character of reconnection dynamics (magnetic structures and 
their motions), and the basic microscopic plasma parameters that determine the radiative 
cooling and the observable radiative signatures. 



-4- 



One issue that complicates the analysis of reconnection dynamics and, at the same time, 
provides a strong reason to believe that reconnection is unavoidable in these systems is the 
expected secondary tearing (aka plasmoid) instability of the equatorial current layer. The 
resulting quasi-p eriodic reconnection and plasmoid ejection have been observed in numer- 
ical s imulations ( iSpitkovskyl l2006t iBucciantini et al.l l2006t iKalapotharakos &: Contopoulos 
20091 ). Indeed, as the equatorial current layer forms just outside the LC and becomes 
thinner and thinner, at some point it will inevitably become unstable to tearing. This 
happens when the aspect ratio of the current sheet — the ratio of it s length to its thick 
ness — exceeds a certain threshold which, as analytical arguments (IBulanov et al.l 11978 



ably of order 100 flSamtanev et al 



Loureiro et al.l 12007; iBhattachariee et al 



2009 



2009 ) and numerical s i mulations suggest, is 



Daughton et al 



2009 



Bhattachariee et al 



prob- 



2009 



Huang fc Bhattachariee 2010; Loureiro et al. 2012; Komissarov et al. 12007; Liu et al. 2011 



Zenitani fc Hoshincl2008l : ljaroschek et al.ll2004l ; Isironi fc Spitkovskvil201ll ; ICerutti et al.ll2012bh 



As a result, the global current sheet breaks into a chain of rapidly growing secondary mag- 
netic islands, separated by much smaller current sheets (see Fig. [2]). Some of these current 
sheets may in turn themselves become so thin and long that they also undergo the same 
tearing instabili ty, forming an hierarchy of smaller and smaller plasmoids and inter-plasmoi d 



current sheets (jShibata fc Tanumal l200ll ; iBhattacharjee et al.l 12009c lUzdensky et al.ll2010l ). 



This hierarchy truncates when the smallest elementary current sheets are no longer unstable 
to the plasmoid instability, i.e., again, when their aspect ratio is about 100. It is these 
elementary current sheets that are the actual sites of magnetic energy dissipation and its 
conversion to plasma energy (and, subsequently, radiation). Determining the basic physical 
plasma parameters in these elementary current sheets is critical for understanding how re- 
connection happens in pulsar magnetospheres and for interpreting its observable high-energy 
radiative signatures. This analysis is the main goal of the present paper. 



This paper is organized as follows. In Section [2] we present the general estimates for the 
physical plasma parameters in the reconnecting CS in the pulsar magnetosphere. In Section [3] 
we apply this model to the Crab pulsar (§ !3.ip and present the physical interpretation of the 
observed high-energy emission components in the framework of the reconnection model. We 
also discuss the conditions for the powerful pulsed synchrotron gamma-ray emission from the 
reconnecting CSs with implications for the observed gamma-ray turn-off line (§ 13. 2p . Finally, 
in Section H] we present our conclusions. 
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pulsar 




reconnecting plasmoid chain 



Light Cylinder 



Fig. 2. — Tearing of the equatorial current sheet in pulsar magnetosphere 



2. Basic plasma parameters inside a radiatively cooled reconnecting current 

layer in the pulsar magnetosphere 



In this section we estimate the basic physical parameters of the plasma inside the smallest 
elementary inter-pl asmoid current layers — those at the bottom o f the reconnection layer 
plasmoid hierarchy (IShibata fc Tanumall2001r lUzdensky et al.ll2010f ). The key distinguishing 
feature of our analysis from traditional reconnection theories is the inclusion of the strong 
optically thin synch r otron r adiative cooli ng of the electrons and positrons inside the layer 
[see also iLyubarskiil (119961 1; iPetril (j2012al )]. An a priori reason why one may expect recon- 
nection in pulsar magnetosphere (e.g., for the Crab) to be in the strong radiative cooling 
regime comes from noting that, for any reasonable plasma density and for a mega-gauss 
reconnecting magnetic field expected at these (light-cylinder) distances, the reconnection 
layer temperature, estimated from the pressure balance across the layer, is so high (in the 
GeV range), that the synchrotron cooling length is very small, measured in centimeters - 
much smaller than the relevant system size. It is thus unavoidable that if reconnection does 
happen in this environment, it must be in the strong-cooling regime. 



We start by outlining the key simplifying assumptions of our model. First, even though the 
plasma is collisionless with respect to Coulomb collisions and hence the particle distribution 
function may be very different from a (relativistic) Maxwellian, for simplicity in this paper 
we will ignore the fine details of the distribution function and will represent the average 
kinetic energy of the particles inside the layer by a single parameter that we will call the 
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temperature T = 7Tm e c 2 , where 7t S> 1 is the thermal Lorentz factor. We are mostly 
interested in the perpendicular (to the upstream reconnecting magnetic field) temperature, 
since this is what enters both the pressure balance across the layer and the synchrotron 
radiative losses (see below). 

Next, we will assume that there is no guide magnetic field inside the reconnection layer. In 
the case of reconnection via an equatorial current sheet of an axisymmetric pulsar considered 
here, this assumption is actually well justified by symmetry considerations — i.e., we are 
dealing here with purely antiparallel reconnection. 

We are interested in reconnection occurring beyond the LC but not too far from it, namely, 
at distances comparable to the LC radius, r — R LC ~ -Rlc- m this region the reconnecting 
magnetic field upstream of the equatorial current layer has both radial and toroidal compo- 
nents which are comparable to each other. In the lab frame, there is also a poloidal electric 
field of comparable strength reversing across the equator. 



The proper way to consider the reconnection problem is to do this in the frame of reference of 
the u pstream plasma, i.e., the plasma just above and just below the current sheet (ILyubarskii 



19961 ). which flies out along rota ting open magne tic field lines with a relatively large bulk 
Lorentz factor of T ~ 10 — 1000 fjLyubarskyl Il995l ) . Indeed, because in the lab (i.e., pulsar) 
frame the perpendicular (to the magnetic field) plasma flow is described by the E x B drift, 
the upstream lab-frame poloidal electric field (which is of order the magnetic field and thus 
is not negligible in the vicinity of the LC) vanishes in this comoving frame. In addition, 
the plasma upstream velocity is, of course, zero in the comoving frame by definition. These 
two simplifications leave us with a clean canonical reconnection problem setup where the 
plasma entering the reconnection layer is static (apart from the flows associated with the 
reconnection process itself). In this approach, the upstream plasma is characterized only by 
three comoving-frame parameters: the upstream reconnecting magnetic field Bq (since there 
is no guide field), the density no, and the pressure P$ (which we shall neglect in our analysis). 
Note that the comoving reconnecting field B should be comparable to the corresponding 
lab- frame magnetic field iLyubarskiil ( 119961 ). This is because, although the bulk Lorentz 
factor of the upstream plasma flow may be large, r > 1, most of it is due to the ultra- 
relativistic parallel motion, which does not change the magnetic field strength, whereas the 
perpendicular plasma velocity, \v±\ = cE po \/B tot = c (r/Rj^c) B po \/ B tot , is only modestly 
relativistic at r ~ R LC (here E po[ and B po \ are the magnitudes of the poloidal components 
of the electric and magnetic field, respectively, and B tot is the total magnetic field strength, 
all in the lab frame). Another way to see that Bq ~ -Bo.iab is to use the fact that B 2 — E 2 is 

On the other hand, we note that 



a Lorentz invariant and hence Bq 



r>2 _ rp2 < r>2 
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the comoving upstream plasma density uq is much lower (by a factor of T) than the density 
in the lab frame. 



Our main goal is to express the physical parameters in the layer — the plasma density n, 
the temperature T, and the thickness of the layer 6, all considered in the comoving frame, - 
in terms of the upstream reconnecting magnetic field Bq, and the reconnection rate, i.e., the 
reconnection electric field E, which are treated h ere a s given known quantities. Our analysis 
is somewhat similar to that of iLyubarskiil ( 119961 ) and iPetril (j2012a ). although there are also 
m any important d i fferen ces (see §H]). It also shares some similarities with the recent work 
bv lArka fc Dubusl boi$ ). 



The main equations governing the physical conditions inside the layer are the following. 
First, in the absence of a guide magnetic field (see above), the comoving frame pressure 
balance across the layer (between the magnetic field outside and the plasma pressure inside, 
since we neglect the upstream plasma pressure) reads: 



d2 

2nT = ^ , 



(1) 



where the factor 2 on the left-hand side accounts for the fact that we have a 2-species plasma 
(electrons and positrons): n e = n e + = n. 



The second equation relating T, n, and 5 is the Ampere's law in the direction of the main 
reconnection layer current (we ignore the displacement current although in principle it may 
lead to a modest correction): 

•OR ° B ° 

j = 2 /3 dr enc ~ — — , 

i.e., 

Bo 



(2) 



(3) 

where cf3 dr is the average drift speed of current-carrying electrons and positrons inside the 
layer; we shall regard (3 d r < 1 as a constant parameter of order unity [in contrast with 
Arka &: Dubusl (120121 ) who assumed it to be small]. We shall also assume that (3 di: is not 
too close to 1, i.e., regard (1 — (3 dr ) to be of order 1 as well. Since the particles are ultra- 
relativistic, /?dr is the magnitude of the average cosine of the angle the particles make with the 
main current direction; it is thus a measure of angular anisotropy of the particle distribution 
function. We note that, as we shall show below, the assumption f3 dT ~ 1 is equivalent to 
the expectation that the layer thickness 5 be comparable to the typical Larmor radius of the 
particles in the layer, which is quite natural and common in collisionless reconnection. 



-8 - 



Finally, the third key ingredient in our analysis is the energy conservation equation in the 
strong cooling regime. This is described as a balance between the electro-magnetic energy 
(Poynting flux) flowing into the layer from both sides, 2 S = 2 (c/47r) EB , and the radiative 



cooling rate per unit area, (see, e.g., Assuming radiative cooling to be 



optically thin, the latter can be estimated in terms of the volumetric cooling rate Q as 
(28) Q ~ (25)2nA(B,T,n), where A(B,T,n) is the cooling function (the radiative power 
per particle). Thus, the overall energy balance reads: 



— EB ~ 2nA(B^ T )8. 

47T 



(4) 



In the case of synchrotron radiation cooling considered here (we neglect IC and curvature 
radiation as cooling mechanisms although they may be important in producing the very high 
energy emission), the cooling function is 



A(5, 7 ) 



2<7 T cf^ 7 2 



(5) 



where B± is the magnetic field component perpendicular to the particle's motion, ot = 
(87r/3)rg is the Thomson cross section and r e = e 2 /m e c 2 ~ 2.8 x 10~ 13 cm is the classical 
electron radius. For a rough estimate, we will assume that most synchrotron radiation is 
radiated by electrons with 7 ~ 7^ in a perpendicular magnetic field Bj_ ~ B [althoug h the 
field may be much lower for the most energetic particles focussed deep into the layer (IKirk 
20041 ; lUzdensky et al.ll201ll ; ICerutti et al.ll2012al )]. Then, our energy balance yields: 

firec 



It ~ 



(6) 



where we introduced a dimensionless reconnection rate parameter 

_ v ICC _ E 

Prec — n i 

C B 

as the ratio of the reconnection inflow velocity, v rec , to the light speed. We shall treat (3 rec as a 
known finite number of order unity. Existing PIC simulations of relativistic pair reconnection 



(7) 



1 Strictly speaking, this is not quite correct since, even in the strong radiative cooling regime considered 
here, a finite fraction of the energy may be carrie d away by the bulk kinetic and thermal energy of the 
reconnection outflow (jUzdenskv fc McKinneyll2011l ): eventually, upon leaving the layer, this energy will be 
converted to the thermal (by viscosity and/or shocks) and residual magnetic energy of the plasmoids; the 
former will eventually be radiated away but because this happens already outside the elementary current 
layer under consideration, this radiation will have no effect on the thermal structure of the layer of interest 
to us here. Thus, the energy leaving the layer results in a modification of our heating/cooling balance by a 
factor of order unity; for simplicity, however, we shall ignore this complication in our present rough analysis. 
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SUffff6St typi-Cc ll Vctlu.6 S Pre.c ^ 

0.1-0.2 ( e.g., IZenitani fc Hoshinol 120071 . 120081 : iLiu et all 1201 ll 
Cerutti et al.ll2012bl ; iKagan et al.ll2012[ ). Note, however, that almost all of these simulations 
did not include the radiation reaction force which, we argue, actually plays a dominant role 
in pulsar magnetosphere reconnection, both as a cooling mechanism and as a strong source 
of effective radiative resistivity. The only study of relati vistic pair reconnection tha t does 
include synchrotron radiation reaction self-consitently, by lJaroschek &: Hoshinol (120091 ). does 
not report a measurement of the reconnection rate. We note, however, that strong plasma 
compression expected in the stron g radiative cooling regime con sidered here may have an 
accelerating effect on reconnection ( Uzdensky fc McKinneylboill ). 



Equations (JTJ) , ([3j) , and (jSJ), represent a system of 3 equations for 3 unknowns n, 7t, and S. 
Solving this system we get: 



It 



n ~ 



ftcc n a 8lTe _oo a PO 

F ~ PrecPdr — O PrecPdr 

GtUQ <JT-D0 r e 

1 Bl 



A/12A- cc /3dr 8irm e c 2 V p ' 



5 ~ — 7tPo ■ 

Pdr 

where we introduced the nominal cyclotron radius: 

,.2 



m P c 



Po 



eB, 



(9) 
(10) 

(11) 



o 



Let us now discuss the physical meaning of these results. 



1 /2 

1) The temperature comes out to be essentially (apart from the factor f3j r of order unity) at 
the classical synchrotron radiation-reaction limit 7 ra d associated with the magnetic field B , 
at which the accelerating electri c force is balanced by radiative losses: eE = eBo(3 rec - 
A(Bn, 7 ra .d)/c = 2a T (BZ/87r) ^ fIGuilbert et al1 ll983l: Ide Jager et allll996l : kharonian etal. 



2002 



Lyutikovll2010l : lUzdensky et al.ll201ll ). Another interesting implication of this estimate 



is that the average projection of the synchrotron radiation reaction force per particle on 
the direction of the reconnection electric field, which can be thought of as the effective 
radiative resistive drag force, is about (3j r times the electric force per particle, eE = eB (3 rec . 
This implies that, even though here we are dealing with a collisionless plasma, the effective 
resistivity due to synchrotron radiation reaction provides an important (of order p^ r ~ 1) 
contribution in balancing the reconnection electric field in the generalized Ohm's law. 



2) As expected in reconnection, apart from factors of order unity, the current layer thickness S 
is basically the relativistic Larmor radius corresponding to this temperature (Q, 5 ~ p(it) = 
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7tPo; this, in turn, equals the total (electrons plus positrons) relativistic collisionless skin- 
depth d = c/ujp = (7t m e c 2 /Rune 2 ) 1 / 2 by virtue of the pressure balance condition (TT]). 



3) The layer's comoving plasma density n is a function of Bq only, independent of the 
upstream density. This is not surprising: in the strong-cooling regime with a small up- 
stream plasma-/?, the upstream ambient plasma entering the layer has to compress by a 
large factor in order to maintain the pressure balance w ith the upstream magnetic field 
Dorman &: Kulsrudl ( 119951 1; lUzdensky &: McKinneyl ( 1201 ll ). In particular, the density inside 
the layer may be significantly higher than the expected upstream comoving density no, in 
which case one expects an enhancement of the reconnectio n rate (lUzdensky fc McKinney 
201ll ). Furthermore, as was suggested by iLyubarskiil ( 119961 ). the upstream plasma density 
may be increased by the pair production due to annihilation of high-energy 7-photons emit- 
ted by the layer. These secondary upstream pairs should be much cooler than the plasma in 
the reconnection layer and could pr ovide a significan t source of optical/UV/X-ray photons 
through their synchrotron emission ( ILyubarskiil 119961 1 . 



It is interesting to note that all three quantities (jt, n, and 5) essentially depend only on the 
strength of the reconnecting magnetic field Bq. It is convenient to normalize this magnetic 
field to a fundamental field scale, namely the classical critical field defined as a field at which 
the nominal po becomes equal to the classical electron radius: 



B, 



m P c 



er P 



2\2 



(m e c 2 ) 



— ~ 6.0 x 10 15 G. 



(12) 



[For reference, this field is a^ 1 ~ 137 times larger than the critical quantum magnetic field 
Bq ED ~ 4.4 x 10 13 G, where «f s = e 2 /hc ~ 1/137 is the fine structure constant.] Then, 
introducing a dimensionless field strength 



b = — - — 

B c \ po 



(13) 



we can rewrite the above expressions as follows (dropping numerical factors of order unity): 



7T ~ (/U9dr) 1/2 &- 1/2 



n 



(/3recA 



drJ 



-1/2 



87rm P c 2 



b 1/2 = (&eM- 1/2 (STrr, 3 )- 1 6 5 / 2 



(14) 
(15) 
(16) 
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3. Observational Implications 

3.1. Application to the Crab pulsar 

Let us now apply these results to the case of the Crab pulsar magnetosphere. We shall 
take the fiducial value of the magnetic field near the LC to be -Boab = 4MG, so that 
6 Crab ~ 6.7 x 10~ 10 and po — r e/b — 4.2 x 10~ 4 cm. We then obtain the following estimates 
for the key reconnection layer parameters (in the wind comoving frame): 

7T,Crab ~ 4x 10 4 (/5 roc /3 dr ) ^ 2 T = lT m e c 2 ~ (/3 rcc /3 dr ) l ' 2 20 GeV , (17) 

ncrab ~ (/3 rec /3 dr )- 1/2 2xl0 13 cm 3 , (18) 
5crab ~ 2 y/fcjK &" 3/2 r e ~ 1.2 x 10 14 r e ^(3 rcc /(3 dr ~ 30 cm V/W&r • (19) 



In the lab (pulsar) frame, the typical energy of the particles should then be in the range of 
hundreds of GeV and the de nsity of order 10 15 cm -3 , assuming the wind bulk Lorentz factor 
in this region of a few tens ( ILyubarskylll995l ). 



Since the reconnection proceeds in the strong synchrotron radiation cooling regime, a sig- 
nificant fraction (~ 1/2) of the magnetic energy entering the reconnecting current sheet is 
promptly radiated away. The overall radiated power should basically scale with the total 
pulsar spin-down power L sd , times the dimensionless reconnection rate parameter /3 rec ~ 0.1 
times some geometrical factors of order unity; overall, it should be as high as a few percent 
of L sd . Since the typical particles in the layer are essentially at the radiation reaction limit 
7 rad , the resulting synchrotron radiation comes out at e p h ~ (3 rP .r 160 MeV ~ 10 — 20 MeV in 



the comoving frame, independent of the magnetic field fide Jager et al.lll996t lAharonian et al. 



20021 ; lLyutikovll2010l ); the emission is then sharply cut off at higher energies. Boosting this by 
a relativistic Doppler factor of a few tens, associated with the overall parallel flow, results in 
an observed radiation in the GeV range. This is in general agreement with the characteristic 
break at a few GeV observed in the pulsed emission of the Crab and a few other pulsars (e.g., 
Abdo et al.ll2010al Jbl). In addition, synchrotron emission in excess of the classical ~ 160 MeV 



limit may be attributed to the efficient focusing of relativistic Speiser orbits of the most ener- 
getic particles deep into the reconnection layer and the resulting suppress ion of synchrotron 



radiation reaction (iKirkl |2004| ; lUzdensky et al.ll201ll ; ICerutti et al.ll2012al ). In any case, we 



feel it is highly plausible that magnetic reconnection in the equatorial current sheet in the 
pulsar magnetosphere (at r ~ Rlc) can be th e main mechanism pow e ring the observed 



(2010 


); 


Arons ( 


2011. 


2012 
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latter assumed that the radiation is produced in the striped wind at much larger distances. 



Furthermore, the typical laboratory-frame energies of the hot (T ~ 10 GeV) electrons and 
positrons in the reconnection current layer are expected to be of order of hundreds of GeV. 
Thus, these particles poss ess energy to produce the very-high- energy (VHE) pulsed emis 



sion recently discovered by I VERITAS Collaboration et al.l ( 120111 ) and MAGIC (lAleksic et al. 



201ll ). This radiation is probably produced by inverse- Compton up -scattering of the UV or 
X-ray emission form the pulsa r, similar to the model proposed by iBogovalov &: Aharonian 
(120001 ); lAharonian et al.l ( 120121 ); however, i n contrast to our mode l, in their picture the wind 
is cold. Alternative views were proposed by lLyutikov et al.l (120121 ) . where the VHE radiation 
is attributed to synchrotron-self- Compto n emission by th e secondary plasma produced in the 
outer gap of the magnetosphere, and by iBednarekl (120121 ) where it is produced by curvature 
radiation. 



Calculating the actual observable light curves of the synchrotron HE emission and the IC 
VHE emission is beyond the scope of this paper. It requires considering the problem in the 
actual three-dimensional geometry of the force-free magnetosphere of an oblique rotating 
pulsar and t aking into account the sp ecial-relativistic light-delay effects, along the lines of 
the study by iBai fc Spitkovskyl (120 lOl ) . 



Finally, it is interesting to think about the implications of the highly-dynamic , non-station ary 
plasmoid-dominated reconnection process for the radio emission of pulsars (lAronsll2012l ). as 
well as for the overall structure of the pulsar wind. First, note that 5 ~ p(7t) ~ 10 cm 
represents the smallest plasma scales in the comoving frame, at the bottom of the expected 
plasmoid hierarchy in the reconnection layer. The corresponding comoving length of the 
smallest current layers is probably about 100 times larger, i.e., on the order of 10 m. In the 
laboratory frame, these scales will be Lorentz-contracted to millimeters and tens of centime- 
ters, respectively. The coherent plasma motions and structures on these scales (and larger, 
since we expect a self-similar hierarchy of plasmoids and current sheets) may contribute to 
powering the radio counterparts to 7-ray pulses (i.e., radio pulses that are in phase with the 
HE pulses), observed in some pulsars, such as the Crab ( lAbdo et al.ll2010bl ). In addition, 
since here we are dealing with the case of reconnection without a guide field, the original 
current layer, as well as the secondary plasmoids (which in three dimensions look like cylin- 
drical flir2^j^p_esJ J _jri^y^ecome disrupted by the relativistic drift-kink instability (RDKI) 



[e.g., iZenitani Sz Hoshinol 120051 ). In principle, one may expect RDKI to cause additional 
dissipation of the residual magnetic energy remaining in the secondary flux ropes and it also 
may introduce additional small-scale substructures contributing to observable radio emis- 
sion. However, although the development and consequences of this instability in the present 
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context of relativistic reconnection with strong radiative cooling are not yet completely un- 
derstood, the only published study on this subject reports tha t tearing instability prevails 
over RDKI in the radiative regime (IJaroschek fc Hoshinoll2009l ). 



As one moves out away from the pulsar, the reconnecting upstream magnetic field becomes 
mostly toroidal and hence the secondary flux ropes and the main reconnection electric field 
become increasingly radially oriented. Furthermore, the secondary plasmoids/flux ropes 
gradually merge (coalesce) with each other, forming bigger structures. Eventually, however, 
in a radially expanding and accelerating relativistic wind, the neighboring ropes lose causal 
contact with each other and the coalescence process stalls, freezes out; then one gets a 
relatively small number of large, mostly radially elongated flux ropes. The resulting geometry 
is favorable for the mechani sm for produc ing the observed bright X-ray knots on the so-called 
Chandra ring, proposed by lAronsI (120121 ). 



3.2. Implications for the Gamma-ray Turn-off Line 



It is interesting to note that not all pulsars detected at radio frequencies also produce pulsed 
(GeV) emission detectable by FERMI. There seems to be a clear boundary in the P — P 
diagram below which radio pulsars do not shine in gamma rays. Following ( jPetril l2012al ) . 
we ask whether one can understand the reasons for this dichotomy in the context of the 
reconnection model. 



As described above, intense gamma-radiation at the characteristic GeV energies (tens of 
MeV in the wind comoving frame) are produced in the reconnection scenario by synchrotron 
radiation from the hot particles in the reconnection layer. This happens when the recon- 
nection process is in the strong cooling regime, i.e., as long as the comoving synchrotron 
cooling time r sync is shorter than the time the particles spend in the reconnection region (in 
the wind comoving frame). Since we expect most of the energy dissipation to take place 
in the inner part of the magnetosphere, within a few light cylinder radii, we can take the 
characteristic length of the reconnection region to be of order /?lc i n the laboratory frame. 
Then, in the wind comoving frame the characteristic time a particle spends in this region 
is about travel ~ -^Lc/cr = P/2tiT. Thus, we can expect a strong HE emission (with a 
luminosity perhaps as high as about 10% of the overall spin-down power) if 

r sync < P/27rr . (20) 



This condition can be recast in terms of the observable pulsar parameters (P, P) as follows. 
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The synchrotron cooling time, using our estimate (Q for the plasma temperature, can be 
evaluated as 

r sync ~7x 10" 9 s (/3 rcc /3 dr r 1/2 , (21) 

where Plc,6 is the characteristic strength of the reconnecting magnetic field near the LC in 
mega-gauss. Comparing with fflUj) , we see that r sync is by a factor (3~l ~ 10 greater than 
5/c. 

Since the magnetic field within the LC is approximately dipolar, B^c can be estimated in 
terms of the equatorial magnetic field at the neutron star surface 5ns as Plc — Pns (Pns/Plc) 3 , 
where Pns — 10 km is the neutron star radius. In turn, Pns is related to the pulsar spin- 
down power L s d = 4tt 2 IPP~ 3 , where / ~ 10 45 gcm 2 is the pulsar moment of inertia, and 
hence to the pulsar spin deceleration rate P as Pns — 2.6 x 10 19 G (PP) 1 / 2 . This results in 
the following estimate of the magnetic field at the LC: 

P LC ~ 3 x 10 8 GP~ 2 ' 5 P 1/2 . (22) 



Substituting this into our expression (J2TJ) for r sync , we can write the condition (120]) for strong 
pulsed GeV gamma-ray emission from the reconnecting equatorial current sheet as 



P > 2 x lO^ 15 T 4 / 3 (/3 rcc /3 dr )- 2 / 3 P 11 / 3 ~ 10" 12 if 3 (& cc (3 di r 2/3 P U/3 
where T2 is T/100, P is given in seconds, and P is dimensionless. 



(23) 



This relationship defines the so-called gamma-ray turn-off line in the P — P diagram. Most 
of the FERMI-detected p ulsars, including m illisecond pulsars, indeed lie on or above this 
line in the P — P diagram (jAbdo et al.ll2010bl ). in agreement with our expectation (123]) . The 
power-law slope in the above relationship (11 /3) is c l ose to what is observed; it is also close 
to the slope <ilogP/<ilogP = 4 obtained by iPetril (j2012a ) [see his equation (14)], which, 
however, appears to be based on somewhat different physical principles. 



4. Conclusions 



In this paper we put forward a basic physical picture of how magnetic reconnection happens 
in the equatorial current sheet just outside of the light cylinder in the pulsar magnetosphere. 
We argued that the current sheet is inevitably unstable to the tearing (aka plasmoid) in- 
stability (and perhaps also the relativistic drift-kink instability), which leads to a violent, 



- 15 - 



highly dynamic magnetic reconnection process involving a hierarchical chain of secondary 
plasmoids/flux ropes of different sizes that are constantly forming and coalescing with each 
other. The actual magnetic energy dissipation takes place in many small elementary inter- 
plasmoid current sheets. Estimating the basic physical conditions within these reconnecting 
current sheets and discussing their observational implications constituted the main objectives 
of this paper. 

We argued that, in strong contrast to reconnection in more familiar solar-system environ- 
ments, relativistically-hot reconnection layers in magnetospheres of pulsars such as the Crab 
are subject to strong radiative cooling. A large part of the magnetic energy flowing into the 
reconnection region as Poynting flux is promptly radiated away through synchrotron radia- 
tion. The resulting heating/ cooling balance condition, in combination with the condition of 
pressure balance across the reconnection layer (in the absence of a guide magnetic field) and 
Ampere's law, enabled us to estimate the key fundamental parameters of the layer — the 
comoving pair plasma density no, temperature T, and the layer thickness 5 — in terms of 
just a single dimensional system parameter, the (co-moving) reconnecting magnetic field B , 
plus a couple dimensionless parameters of order unity, namely the dimensionless reconnection 
rate (3 rec and the average charge-carrier drift velocity in the layer, f3& r . We then proceeded 
to obtain simple scalings of the plasma parameters with the dimensionless magnetic field b 
(i.e., Bq normalized by the critical classical field B c \ ~ 6 x 10 15 G). 

Our estimates show that the comoving plasma temperature inside the reconnection layer is 
essentially comparable to the classical synchrotron radiation-reaction limit (of order 10 GeV 
for the Crab), which implies that, in the wind comoving frame, most of the synchrotron 
radiation is emitted at energies of order (3 rec m e c 2 a" 1 ~ 160 MeV/3 rec , where a ~ 1/137 is 
the fine structure constant. The comoving plasma density inside the layer required by the 
pressure balance is then essentially determined by the upstream magnetic field strength and 
can be much higher than, and relatively insensitive to, the ambient (upstream) density. This 
density enhancement in the layer can be explained, for example, by the strong compression of 
the plasma entering the layer, expected in the strong cooling regime of reconnection. Finally, 
the characteristic comoving thickness 5 of the current layer is expected to be comparable 
with the relativistic Larmor radius of electrons and positrons corresponding to the above 
temperature T (of order 30 cm in the comoving frame in the Crab case). Since the temper- 
ature is close to the radiation- reaction limit, as we discussed above, this layer thickness is 
comparable, essentially by construction, to the corresponding synchrotron cooling length. 

These findings have important implications for our interpretation of the pulsed gamma- 
ray emission observed in many pulsars, including the Crab. First, as was also found by 
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Lyubarskiil ( 119961 ) and IPetril (j2012al ). the main synchrotron emission from the relativistically 
hot reconnection layer in the wind comoving frame is expected to be at around the classical 
synchrotron radiation-reaction limit of e ra( i ~ 160MeV/3rec, where /3 rcc ~ 0.1 — 0.2 is the 
dimensionless reconnection rate. After correcting for relativistic Doppler factor of order 100, 
associated with the overall ambient plasma flow, this synchrotron radiation is boosted into 
the GeV r ange and can thus account for th e pulsed high-energy (GeV) emission detected 
by FERMI iLyubarskiil ( 119961 ); IPetril (j2012al ). Furthermore, by imposing the strong-cooling 
condition (i.e., requiring that in the wind comoving frame the synchrotron cooling time of 
the hot electrons in the layer be shorter than their travel time across the inner part of the 
magnetosphere), we derived the condition [see eq. (|23|) ] for a given pulsa r to be a strong HE 
emitter — i.e., the 7-ray turn-off line in the PP diagram [see also IPetril ( I2012af )]. 



Secondly, the inverse- Compton emission by the same electrons and positrons in the hot 
reconnection layer can produce VHE (hundreds of GeV ) emission di scovered recently i n 
the Crab pulsar by IVERITAS Collaboration et all dioilh and MAGIC JXleksic et al.llioilL 
Finally, the expected collective plasma motions of, e.g., plasmoids/flux ropes in the highly- 
dynamical reconnecting plasmoid chain seem to be in the range of scales at which they may 
contribute to the observed radio emission from the pulsar. 



Our analysis in this paper has a l ot of similarities with those presented by ILyubarskiil (119961 ) 
and, more recently, IPetril ( j2012al ) and lArka fc Dubud ( 120121 ). However, there are also many 
important differences. In particular, we build our model upon the modern understanding of 
reconnection that incorporates the hierarchical picture of reconnection in the plasmoid dom- 
inated regime and the effects of strong radiative cooling and the associated strong plasma 
compression in the reconnection layer. As a result, we are able to obtain explicit approxi- 
mate expressions for the comoving plasma temperature, density, and layer thickness without 
invoking poorly known parameters s uch as the multip li city in the u pstrea m plasma, etc. 



Furth ermore, whereas the studies by ILyubarskiil (119961 ). IPetril (j2012a| ). and lArka Sz Dubus 
( 120121 ) were mostly concerned with explaining the high-energy (GeV) pulsed emission, we 
also addressed the origin of the VHE (hundreds of GeV) and radio emission and argued 
that reconnection in the puls ar ma g netosp here may have important implications for both 
of them. Finally, whereas the IPetril ( 12012al ) model was applied for the pulsar wind at very 
large distances, r 3> Rlc, we develop our model in the context of t he inn er pa rt of the pulsar 
magne tosphere, the near- wind region r ~ -Rlc, as did ILyubarskiil (119961 ) and lArka fc Dubus 



( 120121 ). The authors of the latter paper dArka &: Dubusll2012h . however, are mainly inter 
ested in the case of a non-reconnecting current sheet and do not discuss the reconnection 
rate and the corresponding reconnection-powered plasma heating. Also, they are mostly 
concerned with the case when radiative energy losses from the current layer are small, de- 
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scribing the strong-cooling case (which is the primary focus of our present study) as being 
beyond the scope of their work. Correspondingly, instead of using the energy balance be- 
tween reconnection heating and radiative cooling to obtain a completely closed system of 
equations, they regard the ratio A of the current layer thickness to -Rlc as an arbitrary free 
input parameter (which they typically take to be 0.01). In contrast, in our approach the 
thickness is determined self-consistently by the relevant reconnection-layer physics, including 
the heating/cooling balance. 

We thank Jonathan Arons, Andrei Beloborodov, Benoit Cerutti, Jeremy Goodman, Jason Li, 
and Yuri Lyubarsky for fruitful discussions. This work was supported by NSF grant PHY- 
0903851, DOE Grant Ide^sc0008409l and NASA grants NNX10A039G and NNX12AD01G. 
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